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ABSTRACT 

Aims. The circumstellar dust shells of intermediate initial-mass (~ 1 to 8 M Q ) evolved stars are generated by copious mass loss during the 
asymptotic giant branch phase. The density structure of their circumstellar shell is the direct evidence of mass loss processes, from which we can 
investigate the nature of mass loss. 

Methods. We used the AKARI Infrared Astronomy Satellite and the Spitzer Space Telescope to obtain the surface brightness maps of an evolved 
star R Cas at far-infrared wavelengths, since the temperature of dust decreases as the distance from the star increases and one needs to probe 
dust at lower temperatures, i.e., at longer wavelengths. The observed shell structure and the star's known proper motion suggest that the structure 
represents the interface regions between the dusty wind and the interstellar medium. The deconvolved structures are fitted with the analytic bow 
shock structure to determine the inclination angle of the bow shock cone. 

Results. Our data show that (1 ) the bow shock cone of 1 - 5 x 10~ 5 M Q dust mass is inclined at 68° with respect to the plane of the sky, and (2) the 
dust temperature in the bow shock cone is raised to more than 20 K by collisional shock interaction in addition to the ambient interstellar radiation 
held. By comparison between the apex vector of the bow shock and space motion vector of the star we infer that there is a flow of interstellar 
medium local to R Cas whose flow velocity is at least 55.6 km s -1 , consistent with an environment conducive to dust heating by shock interactions. 



Key words. Circumstellar matter 
kinematics and dynamics 



1. Introduction 



Infrared: stars — Stars: AGB and post-AGB — Stars: individual (R Cas) — Stars: mass loss 



ISM: 



Numerous observations have elucidated the magnitude and ubiq- 
uity of mass loss across the upper-righ t side of the He rtzsprung- 
Russell diagram in the fifty years since lDeutschl (fl956) discussed 
the existence of blue-shifted circumstellar cores in the spectrum 
of the red supergiant star a Her, which constituted one of the first 
pieces of direct evidence for high rates of mass loss. Since such 
high rates of mass loss among asymptotic giant branch (AGB) 
stars rival evolutionary timescales and substantially affect stellar 
evolutionary tracks, careful investigations into mass loss from 
these stars are indeed necessary. Moreover, since mass loss de- 
fines the boundary conditions for stellar evolutionary tracks, the 
rate of mass loss from these stars is hardly time-invariant. These 
facts have had theorists perplexed, who are struggling with the 
basic challenge of how to lift so much mass away from th e grav- 
itational hold of the star (e.g. Gustafsson & Hofner 2004b. 

Observations made with the Infrared Astronomical Satellite 
(IRAS) during the 1980's in the far-infrared (far-IR) have 
demonstrated that extended shells of evolved stars - the an- 
ticipa ted effect of continuous dusty mass l oss - were present 
(e.g. lStencel e7d]|l988tlYoung et al.lll993al) . During the 1990s, 
the Infrared Space Observatory (ISO) and ground-based IR 
work began to refin e those results (e.g. Ilzumiura et al. 1997; 
iMeixner et all [l999), indicating variations in the mass loss 



* Based in part on observations with AKARI, a JAXA project with 
the participation of ESA, and with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of 
Technology under a contract with NASA. 



rate over time that resulted in multiple shells and axisymmet- 
ric structures. This decade, we are fortunate to have new in- 
struments with higher resolution and sensitivity, such as the 
AKARI Infrare d Astronomy Satellite ( AKARI, formerly known 
as ASTRO-F; Murakami et al. 2007) and the Spitzer Space 
Telescope (Spitzer. IWerner et al.1 12004) that can more carefully 
map out the mass loss history of evolved stars. 

In parallel with investigations into the mass loss history of 
evolved stars, evidence of interactions between the circumstel- 
lar matter and interstellar medium (ISM) around AGB stars is 
growing, with ne w observati ons of bow shocks around R Hya 
dlJeta et al.l20"06h and Mira dMartin et alj2007tlUetall2008b, plus 
theoretical considerations of the phenomena (e.g. Villaver et al.l 
2003; Wa reing et al.ll2007l) . This report is one of the first of sev- 
eral studies of well-resolved extended circumstellar dust shells 
of AGB stars under AKARI and Spitzer observing programs la- 
beled "Excavating Mass Loss History in Extended Dust Shells 
of Evolved Stars" (MLHES). In this paper, we explore the de- 
tection of an arcminute-sized dust shell around the oxygen-rich 
AGB star, R Cassiopeiae (HD 224490; hereafter R Cas), in con- 
text of interactions between the AGB wind from R Cas and the 
ISM. 



2. R Cas: the Star and its Circumstellar Dust Shell 

The Mira type variable, R Cas, is an oxygen-ri ch, M7IIIe star 
with a period of 43 1 days dKukarkin et al.ll 19691) and estimate d 
mass loss rate of 5 x 10" 7 M Q yr" 1 dKnapp & Morris 1985). 
This star is known to show an extended circumstellar shell orig- 
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inall y detected by IRAS at 60^m, having angula r extent up to 
43 ([Young et alj|l993al) . lBauer & Stencej d!994l) later reported 
the angular size of 2.'8 using deconvolution of the point-spread- 
function (PSF). The linear extent of the shell depends on the 
distance determination. 

Hipparcos dPerrvman et al.in~997l) measured a parallax for R 
Cas to be 9.37 + 1.10 milli-arcseconds (m as), corresponding to 
107 ± 13 pc. A new calculation done by Ivan Leeuwe 5 (12007ft 
yielded 7.95 + 1.02 mas, corresponding t o 127 + 16 pc. Other 
authors prefer val ues as small as 100 pc dPourbaix et aljr2 003: 
iKnapp et al.l2003b or as large as 220 pc (Knapp & Morris! 1985k 
Hipparcos also detected a 10 mas shift of the centroid of the 
star, whose time-dependent, asymmetric size was determined 
interf erometrically to be 20-40 mas diameter dHofmann et al.l 
1200 lh . Since the rotation of a bright stellar spot, for example, 
could increase the apparent parallax value, the Hipparcos mea- 
surement for the distance t o R Cas is more uncerta in than the 
quoted error values above. IVlemmings et al.l d2003l) performed 
VLBI astrometry and obtained a significantly smaller parallax 
value (5.67 ± 1.95 mas) for R Cas with respect to the Hipparcos 
measurement. Their analysis suggests the distance of 176"'? pc. 
This maser astrometry value is more in line with the larger dis- 
tance esti mate of 160 pc in ferred from the Period-Luminosity 
relations dHaniffetal.1119951) . Since an ensemble of VLBI mea- 
surements over a long time-basis is less affected by uncertainties 
induced by changes on smaller time scales that probably affected 
Hipparcos measurements, we will adopt the VLBI distance mea- 
surement of 176 pc for our purposes. 



3. Far-IR Observations of R Cas 

3.1. AKARI FIS Imaging 

We observed R Cas in the four bands at 65, 90, 140 and 16 0/im 
with the Far-Infrared Surveyor (FIS; Ka wada et alj l2007h on- 
board AKARI on 2007 January 16 as part of the AKARI- 
MLHES Mission Program (PI: I. Yamamura). Observations 
were made with the FIS01 (compact source photometry) slow- 
scan mode, in which two linear strips of forward and backward 
scans were performed with a 70" spacing between the scan strips 
at the 15" s scan rate with a reset rate of 0.5 sec, resulting in 
the sky coverage of roughly 10' x 50' centered at the target. 

The FIS Slow-Scan Toolkit (ver. 20070914*1: IVerdugo et al l 

120071) was used to reduce the data. We found that the quality of 
the resulting map was improved when we used a temporal me- 
dian filter with the width of 200 s (or longer), a temporal boxcar 
filter with the width of 90 s, and the sigma clipping threshold 
of 1.5. For the reduction of data in the short wavelength bands 
(at 65 and 90/im), the results were also improved when we per- 
formed flat-fielding of the data using the local "blank" sky data. 

The resulting maps are in 15" and 30" pixeL 1 (nominal 
scale) for the short wavelength (65 and 90/vm; SW) and long 
wavelength (140 and 160^/m; LW) bands, respectively, with the 
average of 6, 9, 15, and 10 sky coverages per pixel at 65, 90, 140 
and 160/mi, respectively. The resulting 1 <x sensitivities and the 
average sky emission (the component removed by median filter- 
ing during the reduction) are found to be 1.2, 0.6, 1.3 and 1.3 
MJy sr" 1 and 8.7 ± 0.3, 8.6 ± 0.1, 13.2 ± 0.3 and 9.9 + 0.3 MJy 
sr -1 at 65, 90, 140 and 160^/m, respectively. Characteristics of 
the AKARI observations and images are summarized in TableQ] 
images themselves are shown in the top row of Figure Q] 



3.2. Spitzer MIPS Imaging 

R Cas observations at 70/mi by Spitzer was ma de with the 
Multi band Imaging Photometer for Spitzer (MIPS; iRieke et ail 
2004) on 2008 February 18 as part of the Spitzer-MLHES 
project (PI: T. Ueta). Observations were done in the 
photometry/fixed-cluster-offset mode, in which a series of 12 ex- 
posures were made in a spiral pattern around R Cas. This expo- 
sure pattern was intended not to allow the bright central star fall 
on the 70/mi Ge:Ga detector array in "prime" and on the 24fim 
Si:As detector array in "non-prime" in order to avoid saturation 
and severe transient effects due to the central star. In effect, we 
achieved the sky coverage of 13' x 22' centered at the target (but 
the central 3' x 1.8' around the star is unobserved). 

For the data reduction, we started with the basic calibrated 
data (BCD), which are relatively free from instrumental arti- 
facts. However, to optimize the detection of intrinsically faint 
extended shells we used a custom IDL script to remove time- 
and pixel-dependent instrumental effects still remaining in the 
BCD. This script was developed originally to reduce similar 
MIPS data obtained in the MIPS IR Imaging of AGB Dust shells 
(MIRIAD) project (PI: A. K. Speck), based on the idea that it is 
highly unlikely that there arises any periodicity in the time-series 
pixel readings given the way the Spitzer aperture was dithered 
around in the target region. This extra BCD cleaning has proven 
to be effective in correcting those residual instrume ntal effects 
and r emoving residual background sky emission (cf. lUeta et all 

2006) . We then used the MOsaicker an d Point source Extr actor 
(MOPEX) software (ver. 200808193; iMakovoz et al"1l2006l) to 
produce a final mosaicked image. 

The resulting map is in 4'.'92 pixeL 1 (sub-pixelized by a fac- 
tor of 2 from the nominal scale) with the average of 9, and up to 
24, sky coverages per pixel. The 1 <x sensitivity and the average 
sky emission (the component removed during the reduction) are 
found to be 1.1 MJy sr~' and 10.0 + 0.7 MJy sr -1 , respectively. 
Characteristics of the Spitzer observations and image are sum- 
marized in TableQ] image itself is shown in the top-left frame of 
Figure |2] 

3.3. Photometry and Deconvolution 

Photometry was done only with the AKARI images that capture 
the circumstellar shell for its entirety. Because the absolute cali- 
bration of the FIS Slow-Scan data is based on the measurements 
of the diffuse sky emission from zodiacal light and interstel- 
lar cir rus done by CODE/DIRBE measurements ( Verdugo et al. 

2007) , the surface brightness of the extended shell has already 
been calibrated. However, for the emission core that is essen- 
tially a point source a series of corrections (aperture, flux, and 
color corrections) needs to be applied. Thus, we separated the 
structure into two parts - the extended shell and the emission 
core - and performed photometry separately. For the emission 
co re, we followed a met hod of aperture photometry elucidated 
bv lShirahata et al.l d2009l) and evaluated flux correction factors. 
For the extended shell, we simply integrated the surface bright- 
ness over the shell. However, a care was taken not to double- 
count the flux component in the extended shell part that is ac- 
counted for as part of the core emission via aperture correction. 
For each flux value obtained for the core and shell, we indepen- 
dently applied color correction and obtained the final flux values, 
which are listed in TableQ] 



1 Available at http://www.ir.isas.jaxa.jp/AKARI/Observation/ 



2 Available at http://ssc.spitzer.caltech.edu/postbcd/ 
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For PSF calibration purposes, we also ob served an M5 III 
giant B Gru dEngelke. Price. & Kraemerll2006l) with AKARI on 
2006 November 16 and an asteroid Ceres dMiiller & Lagerrosl 
120021) with Spitzer on 2008 February 17. The same instrumen- 
tal/mapping set-up was used for both PSF observations, ex- 
cept for Spitzer observations, for which the photometry/moving- 
cluster-offset mode was used because Ceres is a moving target. 
With the observed PSF maps, deconvolution was performed us- 
ing an IRAff3 task lucy, which is based on the Lucy-Richardson 
algorithm. The average gain and read-out noise were computed 
for each map given the documented detector responsivity for 
FIS and MIPS, integration time, and the averaged number of 
sky-coverage per pixel. Corresponding deconvolved images are 
shown in the bottom row of Figure[T]and in the bottom-left frame 
of Figure |2 



4. Results: the Extended Dust Shell of R Cas 

Both AKARI (Figure Q]) and Spitzer (Figure [2]) images of R 
Cas look very much extended, while the faint surface bright- 
ness in the shell becomes progressively harder to make out in 
the AKARI LW bands, especially at 160yum. The radial surface 
brightness profile centered at the position of the star yielded 3' to 
4' radius for the AKARI SW bands and Spitzer 70/mi band and 
about 2' radius for the AKARI LW bands (at 3cr sky ; Table 
While the position of the star is obviously off-centered, the ex- 
tended shell is roughly circular whose 3<x sky radius is 2'.3 to 2.'8 
at 65, 70, and 90//m. To clarify we also note that the AKARI SW 
band images (N60 and WIDE-S) are affected by the cross-talk of 
the FIS detector that manifests itself as a linearly extended emis- 
sion structure emanating from t he central source into the position 
angles at 48° and 228° E of N (IShirahata et alj|2009b . 

The shell's emission structure consists of the relatively flat 
"plateau" region on the west side (of surface brightness ~ 15 to 
20 MJy sr -1 at 65 to 90//m and < 10 MJy sr~' at 140jum) and 
the region of higher surface brightness (emission core) on the 
east side around the central star. This particular emission struc- 
ture, however, does not appear to be caused by the off-centered 
central star. The deconvolved maps in the AKARI SW bands 
(the bottom-left images in Figure Q]) and in the PSF-subtracted 
Spitzer JOfim band (the bottom-left image in Figure [2]) all show 
consistently that the surface brightness is enhanced along a rel- 
atively well-defined rim that goes around the periphery of the 
shell for almost half the shell on its east side. 

Assuming that the detected far-IR emission is mostly of dust 
continuum, the surface brightness enhancement on the east side 
can be caused by either the density or temperature enhancement 
of dust grains in the shell (or both). This is because in such an 
optically thin environment at far-IR the surface brightness S v is 
proportional to t v B v (T), where r v is optical depth along the line 
of sight and B y (T) is the blackbody function for dust grains at 
the temperature T. Thus, we fit the t v B v (T) curve with the mea- 
sured surface brightnesses at 65, 70, and 90/im to infer the dust 
temperature T^ust and the optical depth at 70yum, tjq^ simultane- 
ously at each pixel. This was done by rescaling the Spitzer 70/jm 
map to the same pixel scale as the AKARI SW band images (15" 
pixeL 1 ) and performing the fitting using PSF-subtracted maps. 
The results are shown in the right frames of Figure[2] 



3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



The dust temperature map (the top-right frame of Figure [2]) 
shows the peripheral enhancement similar to the one seen in 
the (deconvolved) surface brightness maps (the bottom frames 
of Figure [TJ. The dust temperature is the highest (> 20K) in 
the well-defined peripheral rim on the east side of the shell and 
in the less well-defined, clumpier structures on the west side. 
However, the map (the bottom-right frame of Figure [2]) 
does not show any obvious enhancement near the periphery of 
the shell. This suggests that the observed surface brightness en- 
hancement in the extended shell of R Cas is due to the temper- 
ature enhancement in the shell rather than the density enhance- 
ment. What would cause such a temperature enhancement in the 
she ll, then? 

IVlemmings et al.l ((2003) measured proper motion of R Cas, 
<jj. a ,Hs) = (80.52 + 2.35, 17.10 ± 1.75) mas yr" 1 , by VLBI as- 
trometry. This translates to the position angle of 78°0 ± 1°6 east 
of north, which agrees remarkably well with the direction along 
which there is a positive gradient of surface brightness. In other 
words, the emission structure of the extended shell appears to 
show mirror symmetry with respect to a line defined by the di- 
rection of proper motion. Interestingly, the direction of the mea- 
sured proper motion is parallel to the direction of the apparent 
shift of the central star from the geometric center of the shell. 
The shape of the circumstellar shell was fit by an ellipse using 
the MIPS 70pm map, in which the pixel resolution is the high- 
est. We defined the shell edge at the 3cr sky level and measured 
the radius (the distance from the star to the edge) in all direc- 
tions at a certain azimuthal interval. The best-fit semi-major axis 
length (a) and eccentricity (e) were thus searched for by adjust- 
ing these two values while minimizing the difference between 
the measured shape and the ellipse defined by the particular (a, e) 
pair. During this fit, we assumed that the central star would fall 
on one of the foci and the semi-major axis would lie along the 
vector defined by the direction of the measured proper motion. 

The best-fit semi-major axis and eccentricity pair, (a, e), 
turned out to be a = 165V3 and e = 0.3 (hence, the semi- 
minor axis length b = 157'.'7): indeed, the extended shell of R 
Cas is not quite circular in projection. According to this best- 
fit ellipse, the distance from the ellipse center to one of the 
foci is c — a x e — 48'.'8. This means that the central star is 
48" 8 displaced from the ellipse center over the course of the 
shell expansion. At the preferred distance for R Cas, 176 pc 
(IVlem mings et al. 2003), the semi-major and semi-minor axes 
correspond to 0.13 to 0.14 pc. If we assume that the measured 
CO e xpansion rate of 12 km s _1 (Buiarrab al. Fuente. & Omontl 
1994) is the (constant) expansion velocity of the extended dust 
shell, the crossing time of the shell is therefore roughly 10 4 
years. 

If the elongation of the R Cas shell is solely due to the mo- 
tion of the central star with respect to the shell (in an otherwise 
stationary local environment), the star must have been moving 
roughly at 5 mas yr -1 , which is much less than the observed 
proper motion of 82.3 mas yr -1 . Thus, the elliptical elongation 
of the shell and the offsetting of the central star do not seem to 
be self-inflicted, as in the case where a pile-up of AGB wind ma- 
terial occurs at the interface between fast a nd slow AGB wind s 
and defines the edge of the observed shell (Ste ffen et al.| [l998). 
Rather, we speculate that the shell shaping is orchestrated by the 
interactions between the AGB winds and ambient ISM, as in the 
AGB wind-ISM interaction disco vered around an AGB star, R 
Hya, by Spitzer dUeta et al.ll2006l) . In this scenario, the temper- 
ature enhancement seen in the east side of the shell is due to 
dust heating along the contact discontinuity between the AGB 
wind and the ISM flow. The temperature enhancement seen in 
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the west side is then probably due to dust heating in the wake 
of the AGB wind-I SM interactions flowing downstream (e.g. 
IWareing et al.l2 007). This dust heating is likely due to collisions 
rather than shock emission because t here is no known Ha emis- 
sion source at the position of R Cas (Finkbeiner 2003). This in- 
terpretation is consistent with recent H I observations of R Cas 
( Matt hews & R eid 2007), which identifies a head-and-tail struc- 
ture whose alignment in the western wake of our far-IR maps is 
nearly perfect . 

Unlike the R Hya case, however, we do not see any clear 
parabolic bow shock structure that is expected to arise at the 
interface between the AGB wind and ISM. Still, we recognize 
rather round temperature-enhanced regions toward the windward 
direction in the shell and wake-like temperature-enhanced struc- 
tures toward the leeward direction. This is similar to the bow 
shock of Betelgeuse, for which its bow shock appears rather cir- 
cular due to the inclination angle of the shock surface (lUeta et al.l 
2008). Assuming momentum conservation across a physically 
thin (i.e. radiative-cooling dominating) shock interface between 
the AGB wind and ISM, one can express the bow shock shape 
analytically as a function of the latitudinal angle 9 measured 
from the apex of the bow with re spect to the position of the cen- 
tral star as follows (Wil kinll 19961) : 



R(0) = R Q 



V3(l -6>cot6>) 
sin 9 



(1) 



where Rq is the stand-off distance between the star and bow apex 
defined as 



Rn = 



(2) 



for which M is the rate of mass loss, v w is the isotropic stel- 
lar wind velocity, pism is the ambient ISM mass density, and v» 
is the space velocity of the star. Because the observed enhance- 
ment in the shell (i.e. an inclined bow shock) represents a conic 
section of the bow shock paraboloid (in optically thin dust dis- 
tribution the column density becomes the greatest at the cross 
section of the structure with the plane of the sky), we can de- 
termine the inclination angle by fitting the apparent shape of the 
bow s hock with the an alytical function and a given inclination 
angle dUeta et al.ll20"08h . 

The best-fit of the Wilkin solution fitting yields the inclina- 
tion of the bow shock to be 6>i nc i = ±(68° + 2°) and the posi- 
tion angle of 74° + 2°. The inclination angle has a degeneracy 
because this fitting method alone would not determine if the 
structure is inclined to us or away from us with respect to the 
plane of the sky. Nevertheless, the direction of the relative mo- 
tion of the AGB wind shell with respect to the ambient ISM 
determined from the apparent bow shock structure is very much 
different than inferred from the observed motion of the central 
star (6*; nc i = 18? 4 toward us at the position angle of 78°0). The 
discordance between these two values can be resolved only by 
incorporating the idea that the ISM itself is flowing in a partic- 
ular direction. In other words, the difference between the vector 
of the apparent (heliocentric) space motion of the central star 
v* jQ and the vector of the space motion of the star relative to the 
ambient ISM v.jsm (as derived from the 3-D orientation of the 
shock structure) is the heliocentric flow vector of the ambient 
ISM Vism, q : 

VlSM,© = y*,Q - V»,ISM- (3) 

The stand-off distance of 1.'4 ± O.'l translates to O.lpc at the 
adopted VLBI-measured distance of 176pc. Solving eq.® for 



v*,ism with known quantities (M = 5 x 10~ 7 M yr , v w = 12 
km s _1 , and pism ~ l-4mH«isM' where 1.4m H is the average ISM 
particle mass and «ism is the ISM number density), we have 



VlSM.O - 



"ism,© 

^ISM,0 
V ^ISM.O > 



-22.9 ^ 
+67.4 
+ 14.3 



1 



V n ISM 



±33.0 ^ 
+ 12.6 
+3.6 



(kms- 1 ) (4) 



as a function of the remaining unknown, «ism- The Vtsm,o-"ism 
relation suggests that vism,© reaches its minimum of 55.6 km s _1 
at «ism = 0.58 cirT 3 when the bow shock cone points to us and 
72.5 km s _1 at «ism = 75.5 crrT 3 when the bow shock cone 
points away from us. In general, high values are unlikely for 
both Vism,© and «ism- Therefore, the case where the bow shock 
cone points away from us seems unlikely. If this is true the di- 
rection of the ISM flow should be going into the plane of the sky 
(positive v ra d), given the space motion of the central star. Hence, 
the ambient ISM number density needs to be less than 2 cm -3 , 
which is in line with the fact that a relatively higher «ism value 

is expected for R Cas, a low galactic latitude object; 

In the frame of stellar winds at 12 km s _1 (Knapr Tet al.l 

1998), the relative velocity of the ambient ISM flow into the 
interface regions is 12 + 35.5/ V"ism k m s -1 , which translates 
into at least 37 km s if «ism > 2 cm -3 . Such a velocity could 
induce relatively weak (collisional) shock interactions that can 
raise the dust temperature of the AGB wind-ISM interface re- 
gions to ~ 20 - 30 K as seen from the present far-IR data. Given 
the nature of the central star, it is highly unlikely that dust heat- 
ing is solely due to radiation from the central star. Typically 
the ambient interstellar radiation field is thought to contribute 
in heating dust grains in the outermost regions of the extended 
circumstellar shells. The present analysis of the far-IR data in- 
dicates that the AGB wind-ISM interactions can contribute as 
yet another source for dust heating in an otherwise cold envi- 
ronment. In future radiative transfer studies with our far-IR data 
sets, we will quantify relative contributions of these dust heat- 
ing sources at the periphery of the extended dust shells around 
evolved stars as they would impact the physical and chemical 
conditions of the AGB wind ejecta upon injection into the ISM. 

The dust temperature and optical depth maps together with a 
surface brightness map would yield crude estimates of dust mass 
in the shell: 



2r>2 



M dust 



F V A Z D 



(5) 



where F v is the measured flux at a given frequency v of the 
band, A is the wavelength of the band, D, is the distance to the 
source, k is Boltzmann constant, k v is the dust opacity, A\ is the 
area subtended by each pixel, and the summation over pixel i 
refers to integration over the shell. Here, we assume a spheri- 
cal gr ain (0.1/im radius) of olivine (MgFeSiCU; Dorschne r et al.l 
1995) and calculated its opacity with Mie theory (48, 41, and 25 
cm 2 g -1 at 65, 70, and 90/im, respectively) and use the adopted 
distance to R Cas of 176 pc. Our data then suggest that the dust 
mass in the shell is roughly 1 - 5 x 10~ 5 M Q . Given the mass 
loss rate of 5 x 10~ 7 M yr -1 and the dynamical age of the shell 
of ~ 10 4 years, the rough estimate of the total mass in the shell 
is 5 x 10~ 3 M Q , yielding the gas-to-dust ratio between 100 and 
500. The amount of total mass estimated this way is a lower 
limit because the duration of mass loss could have been more 
than 10 4 years. Then, the estimated gas-to-dust ratio will be cor- 
respondingly higher, indicating that there exists less dust in a 
colder environment of the extended dust shells. However, this is 
based on a very rough estimate and needs to be confirmed by 
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radiative transfer calculation incorporating at least dust heating 
by the interstellar radiation field. 

5. Conclusion 

We obtained far-IR images of an oxygen-rich Mira variable R 
Cas at 65, 70, 90, 140, and 160/im using AKAjRJ and Spitzer and 
revealed its very extended (2' to 3' radius, corresponding to 0.1 
pc at its adopted distance of 176 pc), slightly elliptical (e = 0.3) 
dust shell, in which the central star is located offset from the geo- 
metric center of the shell in the direction of the measured proper 
motion of the central star. We recognize a positive gradient of 
the surface brightness along the direction of the "shift" of the 
central star, which apparently is caused by the surface bright- 
ness enhancement along the periphery of the shell in the east 
side seen in deconvolved images. 

Fitting of the surface brightness using data in the 3 shortest 
bands suggests that the observed enhancement is caused by the 
temperature enhancement rather than the density enhancement, 
prompting a need to warm up dust grains primarily on the east 
side of the outer rim. Given the coincidence between the direc- 
tion of the proper motion of the central star and the direction 
of the apex of the peripheric surface brightness/temperature en- 
hancement in the shell, we infer that the observed shell structure 
represents the contact surface of the AGB wind-ISM interaction 
which is inclined to give an overall spherical shape instead of 
a typical parabolic bow shock structure. This AGB wind-ISM 
(collisional) interaction therefore warm up dust grains in the in- 
terface regions, causing the temperature enhancement towards 
the windward direction of the shell. Using the maps, we also 
estimated the total dust mass in the shell to be 1 - 5 x 10~ 5 M . 

The shape of the observed enhancement was fitted with the 
analytical function for the bow shock cone to derive the incli- 
nation angle of 68°. The apex vector of the bow shock cone and 
the space motion vector of proper motion of the central star were 
compared to deduce that there is an ISM flow local to R Cas 
that has a flow velocity of at least 55.6 km s _1 . Then, the rela- 
tive velocity of the ambient ISM flow with respect to the AGB 
wind-ISM interface regions is at least 37 km s _1 . Therefore, such 
weak shocks can play a role in heating dust grains in the outer- 
most regions of these extended dust shells around evolved stars 
in addition to the interstellar radiation field that is originally ex- 
pected to play a role in an environment where luminosity from 
the central source is not enough for required dust heating. 
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Table 1. Characteristics of the Far-IR Observations of R Cas with AKARI and Spitzer 



Date & Time 




A 




Pixel Scale 


Sky Coverage 


S i/ ; sky 


C"sky 






(UT) 


Band 


(jim) 


(jum) 


(arcsec) 


(pixeP 1 ) 


(MJy sr" 1 ) 


(MJy sr" 1 ) 


(arcmin) 


(Jy) 


AKARI/FIS Observations 


2007 Jan 16 10:19 


N60 


65 


22 


15 


6 


8.7 


1.2 


3.1 


121.48+19.90 




WIDE-S 


90 


38 


15 


9 


8.6 


0.6 


3.8 


73.42±10.03 




WIDE-L 


140 


52 


30 


15 


13.2 


1.3 


2.0 


25.89±2.09 




N160 


160 


34 


30 


10 


9.9 


1.3 


1.8 


22.41±7.47 


Spifzer/MIPS Observations 


2008 Feb 18 04:02 


MIPS70 


70 


19 


4.92 


9 


10.0 


1.1 


3.5 






RA offset (arcmin) 

Fig. 1. AKARI/FIS maps of R Cas in the SW bands - N60 (65//m) and WIDE-S (90fxm) at ^"pixeF 1 scale - and in the LW bands - 
WIDE-L (140//m) and N160 (160yum) at SCTpixeF 1 scale - in the top row from left to right, respectively. Background emission has 
been subtracted by a combination of temporal filters during data reduction. RA and DEC offsets (with respect to the stellar peak) 
are given in arcminutes. The number at the bottom left indicates the peak surface brightness in MJy sr -1 . The log-scaled grayscale 
surface brightness is contoured at 90, 70, 50 (in white), 30, 10, 3, and lcr^y. North is up, and east to the left. Images in the bottom 
row show deconvolved maps at corresponding wavelengths, for which "LR" on the bottom right indicates the Lucy-Richardson 
algorithm used. 
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2 0-2-4 2 0-2-4 

RA offset (arcmin) 



Fig. 2. Spitzer/MIPS map of R Cas in the 70yum band at 4"92 
pixel -1 scale (top left frame) and its PSF-subtracted image (bot- 
tom left frame). Background emission has been subtracted by a 
combination of temporal filters during data reduction. RA and 
DEC offsets (with respect to the position of the star, indicated by 
the "star" symbol) are given in arcseconds. Also shown are the 
dust temperature (top right) and optical depth (at 70/mi; bottom 
right) maps derived from the fitting of the shortest 3 waveband 
images (65, 70, and 90//m) under the assumption of optically 
thin dust emission. The wedges on the right indicate the log- 
scaled dust temperature in K (top) and optical depth at 70/im 
(bottom). North is up, and east to the left. 
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